o 



CO 
CN 



Spin waves in diluted magnetic quantum wells 

P. M. Shmakov 2 , A. P. Dmitriev 1 ' 2 , and V. Yu. Kachorovskii 1 ' 2 

1 Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 

2 A.F.Ioffe Physico-Technical Institute, 194021 St. Petersburg, Russia 

(Dated: July 27, 2010) 

We study collective spin excitations in two-dimensional diluted magnetic semiconductors, placed 
into external magnetic field. Two coupled modes of the spin waves (the electron and ion modes) 
are found to exist in the system along with a number of the ion spin excitations decoupled from 
the electron system. We calculate analytically the spectrum of the waves taking into account the 
exchange interaction of itinerant electrons both with each other and with electrons localized on 
the magnetic ions. The interplay of these interactions leads to a number of intriguing phenomena 
including tunable anticrossing of the modes and a field-induced change in a sign of the group velocity 
of the ion mode. 
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Diluted magnetic semiconductors (DMS) have 
recently been the subject of great interest [l|, [2[ 
due to their potential in combining magnetic 
and semiconductor properties in a single mate- 
rial. The DMS are formed by replacing of cations 
in ordinary semiconductors with magnetic ions, 
typically Mn ions. Strong exchange interaction 
between the itinerant electrons and the electrons 
localized on d-shells of the magnetic ions leads 
to a number of remarkable features of the DMS. 
In particular, it results in the effective indirect 
interaction between the ion spins thus promis- 
ing for creating room-temperature ferromagnetic 
systems that may offer advantages of semicon- 
ductors. It also dramatically enhances the effec- 
tive coupling of the itinerant electrons with the 
external magnetic field. In contrast to conven- 
tional GaAs/GaAlAs systems, where small val- 
ues of g— factor prevents manipulation of the spin 
degree of freedom, the giant electron Zeeman 
splitting arising in the DMS as a manifestation 
of the exchange interaction can be on the order 
of the Fermi energy [3|, |j| , offering a wide range 
of spintronics applications. 

Here we discuss spin excitations in the two- 
dimensional DMS. Our studies are motivated by 
recent experiments [5|— l7|] and a theoretical dis- 
cussion [7H9( focused on the spin dynamics in 
diluted magnetic Cdi-zMn^Te quantum wells 
placed into the magnetic field [lj. In Ref. [5J, 
the spectrum of the spin waves, w(/c), was mea- 
sured. Only one excitation mode was observed. 
It was demonstrated that the excitations exist in 
a finite range of wavelengths, k < fc m , and their 
group velocity is negative: dco(k)/dk < 0. The 
experimental data were interpreted [5j,l8| in terms 



of conventional spin waves in the Fermi liquids 
11| , while k m was attributed to the edge of the 



Stoner continuum of the single-particle spin exci- 
tations. Such interpretation implies that the only 
effect of the magnetic ions on the electron spin 
waves is the strong renormalization of the elec- 
tron Zeeman splitting. However, more recent ex- 
perimental observations [g, [7[ supported by the- 
oretical studies 0, |9j appear to be in disagree- 
ment with this conclusion. Indeed, in Refs. [fj, 0] 
two modes of the collective homogeneous (fc = 0) 
spin excitations were observed in Cdi-^Mn^Te 
wells. The modes were identified [g, [D, |9j as 
the spin excitations of delocalized electrons (the 
electron mode) and the electrons on d-shells of 
Mn ions (the ion mode). The dependencies of 
the frequencies 0^2(0) of observed modes on B 
are shown schematically in Fig. ([T]). The most 
important observation is the anticrossing of the 
modes which occurs at a certain "resonant" value 
of magnetic field B = B res [fj, |7| . As it was also 
shown in Ref. J7J , other types of the homogeneous 
spin modes may exist in the system correspond- 
ing to excitations of the ion spins decoupled from 
the spins of the itinerant electrons. 

In this paper, we develop a theory of the spin 
waves in diluted magnetic quantum wells placed 
into magnetic field. We study analytically two 
collective modes which correspond to coupled 
propagation of the electron and ion spin excita- 
tions. We also discuss the ion modes decoupled 
from the electron system. To describe the homo- 
geneous spin oscillations (k = 0), it is sufficient to 
take into account only one type of exchange inter- 
action: the interaction of the itinerant electrons 
with electrons localized on the Mn ions. The thus 



obtained results coincide with those presented in 
Refs. [6J, lZ|. For k ^ 0, the exchange interaction 
between delocalized electrons comes into play. 
Our main purpose is to demonstrate that the si- 
multaneous presence of two types of exchange in- 
teraction gives rise to interesting phenomena the 
most remarkable one is the magnetic-field-driven 
anticrossing of the spin modes. In contrast to 
the case k = 0, the anticrossing can take place in 
a wide range of the magnetic fields and may be 
tuned by the field to occur at a certain value of 
k (see Fig. 2). 




FIG. 1: Anticrossing of the electron and ion spin pre- 
cession frequencies Q e and Qj at "resonant" mag- 
netic field B — B reB ; fi e is the edge of the Stoner 
continuum. 



We consider the 2D degenerate electron gas in- 
teracting with randomly placed magnetic ions. 
The electrons are located in the r = (x, y) plane 
and occupy the lowest level in the well. The 
ions are distributed homogeneously with the 2D 
concentration nj, which is assumed to be much 
higher than the electron concentration n e . The 
magnetic field is applied parallel to the well plane 
(B || e x ). The field leads to the Zeeman split- 
ting of the electron and ion spin levels with en- 
ergies hLj e and HQj, respectively, while the or- 
bital motion remains intact. The Hamiltonian 
of electron-ion exchange interaction reads Hj e = 
-a/2j2k^^kS(r - r fc )|*(zfc)| 2 , where <x is the 
Pauli matrix vector, J^ are the spin operators 
of the ions located at the points R k = {r k ,Zk), 
and &(z) — y /, 2/asm(irz/a) is a wave func- 
tion of the lowest level in a rectangular well 
of width a. Since nj 3> n e , the distance be- 
tween the ions is much smaller than the electron 
wave length and the mean-field approximation is 
applicable. In this approximation, we first re- 
place <xJfc with (o - ) Jfe + o"(Jfe), where averag- 
ing is taken over density matrix of the system 
p. After such decoupling one may search the so- 



lution of the quantum Liouville equation for p as 
a product of the electron and ion density matri- 
ces: p — p e (r, r', t) Y\ k Pk{t)- The averaged spins 
of the electrons and ions are given by s (r, t) — 
n- 1 Js(r,p,t)d 2 p/(27Thl 2 , 3 k (t) = Tr(p k (t)3 k ), 
where s(r, p,i) = Tr(<x/)/2 is the electron spin 
density and / = /(r, p, t) is the Wigner function 
corresponding to p e . As a next step, we replace 
3 k (t) with a smooth function J(r, z, t). Doing so, 
one finds the electron spin precession frequency 
in the ion-induced exchange field 



u> e j(r, t) — cm ,/J(r, z, t)/ha, 



(1) 



[here J(r, t) = J dz\^(z)\ 2 3(r, z,t)\, and the local 
frequency of the ion spin precession 

u Je (r,z,t) =an e s Q (r,t)\y(z)\ 2 /h. (2) 

In addition to the electron-ion exchange inter- 
action, we take into account the isotropic ferro- 
magnetic electron-electron exchange interaction 
by adding the following term 11 | 



u ee (r,t) = -2Gn e v 1 s (r,i), 



(3) 



to the electron spin precession frequency. Here 
G < is the interaction constant (we assume 
\G\ < 1), v = m/2irfi 2 , and m is the electron 
effective mass. 

The equilibrium electron spin den- 



sity, 



^<i( 



[n t (e) - raj. (e)] e x /2 (he 



n u (e) = {exp[(e T M}°/2 - E F )/T] + 1} 1 and 
e = p 2 /2m), is expressed via the effective Zeeman 
splitting fi-Slg. The frequency Q® = uj e +u> e j+uj ee 
is found self-consistently from Eqs. (flj and ([3]): 

n2 = n e /(i + G), (4) 

where 



O e = oj e + anjj° q /ah, 



(5) 



is the effective electron Zeeman splitting renor- 
malized by exchange interaction with the ions. 
In deriving these equations, J(r, z) was sub- 
stituted with the equilibrium ion polarization, 
J e i = jpe x = lB 5/2 {hrij/T)e x , where Bj(x) is 
the Brillouin function. We also assumed that the 
equilibrium exchange field acting on the ions is 
small, an e s x q /a sw ahQ®n e /4EFa ^C hflj, which 
implies that the equilibrium ion polarization is 
not affected by exchange interaction. In con- 
trast, the electron Zeeman splitting is strongly 



enhanced due to high ion concentration, so that 
il e ^$> | w c | (w e < 0, because of the negative elec- 
tron g-i actor, which explains non-monotonic de- 
pendence of n, e on B shown in Fig. 1 [a, lZ|). 

The out-of-equilibrium spin dynamics can be 
described by the Landau-Silin equation [12| for 



the electrons and the Bloch equation for the ions: 



dt 



^'-^IH ly1 =°' w 



[ftje x + U),Je] X J = 0. 



(7) 



Here [■■•] and {•■■} stand for the commuta- 
tor and the anticommutator, respectively, and 
£ = — h[uj e e x + uj e j + cj ee ]&/2. For Ml° e <C Ep, 
Eqs. ([5]) and ([7]) give the following system of 
equations for the perpendicular (with respect to 
B) components of the electron and ion spins 

(js — — 

— + (> F nV + ifl° e )(a + Gs ) = 6 1 (iJ + £nVJ), 



dJ -o 7 

— +injj 

at 



iS 2 s . 



(8) 



Here s = s y + is z , J 



J y + iJ z , vf is Fermi 



velocity, n = (cos ip>, simp), ip is the velocity an- 
gle in the well plane, 6\ — anjs x q /ha, 6 2 — 
3an e J x q /2ha, £ — vf/^1 and s — f Q sd(p/2TT 
[13| . After Fourier transform of Eqs. (jHJ, one 
can find the dispersion equation for the collec- 
tive modes 



'1- 



u 6 2 + G9? e (uj - flj) 

(uj - n° e ) 2 ~ uj-n° e 6 2 + n e (uj-rij) 



vpk 2 



where 6 = \ / 6i6 2 . For k = 0, we get @, Q, 

« ll2 (o) = (n e + Qj)/2 ± ^/(n e -nj) 2 /4 + 6 2 . 

The anticrossing occurs when Q e (B) = Qj(B). 
We see that the constant G drops out from 
011,2(0). In contrast, the dispersion of the collec- 
tive modes strongly depends on the relation be- 
tween \G\ and the dimensionless parameter 8/Q e . 
For 6/fl e <C \G\ (this was the case in the experi- 
ment [a]), the anti-crossing occurs for B < B res , 
when il e > ilj (see Fig. 2a). To see this, one 
may consider the case 6 = (coupling between 
transverse components of electron and ion spins 
is turned off) as a first approximation. This ap- 
proximation was implicitly used in Ref. [8|. For 
5 = 0, there are two branches of the spectrum 
(dashed lines in Fig. 2a), corresponding respec- 
tively to the Fermi-liquid spin waves with nega- 
tive dispersion and the dispersionless excitations 
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FIG. 2: (a) Anticrossing of the two collective spin 
modes for B < B re3 ; (b) Sign inversion of the group 
velocity of the ion mode for B > B. 



of the ion spins. Importantly, for B < B res these 
two branches intersect each other. Turning on 
a finite coupling between modes, 5^0, results 
in the anticrossing, which, for B close to B res , 
occurs at the point 



k res = y/2\G\il e (il e -nj)/v F (l + G). (9) 

Remarkably, k res depends on B, so that the an- 
ticrossing position may be tuned by the exter- 
nal field. The splitting between the modes for 
k Kt k res is given by uji(k res ) - uj 2 (k res ) w 26. 

As seen from Fig. 2a the upper branch of the 
spectrum at a certain wave vector k m reaches the 
Stoner continuum (single-particle excitations), 
which is defined by inequality \u — £l®\ < VfH. 
For k > k m the corresponding ion-type excita- 
tions slowly decay in time due to weak exchange 
coupling with the system of itinerant electrons. 
This decay is similar to the well-known Landau 
damping in plasma [12l j , so that the decay rate 7 
is calculated in a quite analogous way, yielding 



7; 



6 2 n J v F y / k 2 -k 2 



fi°[(l + G) 2 v 2 F (k 2 - k 2 n ) + G 2 n 2 



(10) 



As a function of k, 7 has a maximum. The 
maximal value is given by j max = 6 2 /2\G\il e ~ 
3an e s x q /Ah\G\a. Using data of Ref. 7] (n e = 
0.7 • 10 n cm- 2 , a = 80A, a = 1, 5 • 10" 23 eVcm 3 , 
\G\ « 0.2 s% q w 0.2) we find 7max « 2 • 10 9 s- 1 . 

Another interesting phenomenon arising due 
to simultaneous presence of two types of inter- 
action is a change in a sign of the group ve- 
locity of the ion mode. It can by understood 
by analyzing the spectrum in the limit k — ¥ 0, 
when uo 12 (k) ~ Wi 2 (0) + v 2 F k 2 /2/3. Here /3 = 

n°M, 2 (o)-n e }[(uiAQ)~n e ) 2 +5 2 }/uj 1 , 2 (o)6 2 ,so 

that the dispersions of the modes are controlled 



by signs of Wi(0) — 0° and W2(0) — fig, respec- 
tively. As seen from Fig 1 , there is a critical filed 
B = B , at which Wi(0) = Q°. For B < B Q , 
both spin-wave branches are below the Stoner 
continuum and have negative dispersions. While 
B increases, approaching Bq, the ion spin- wave 
branch becomes shorter (k m — > 0) and disappears 
when B = Bq. For B > Bq, this branch appears 
again above the Stoner continuum and has a pos- 
itive dispersion as shown in Fig. 2b. The disper- 
sion of one of the modes can also change sign for 
Q e > Q,j, provided that <5/O e 3> \G\. 

Above we discussed the coupled collective 
modes. Now we notice that Eqs. ([8]) have also 
a solution s = and J ~ 0. Importantly, the 
latter equation apart from the trivial solution 
J(r, z,t) = has also non-zero solutions obey- 
ing the constraint J dz\ 1 $(z)\ 2 J(r, z, t) = 0. Such 
solutions were called "decoupled" modes [7]. To 
find a number of such modes one should go be- 
yond continuous approximation and replace inte- 
gration over dz in all above equations with the 
sum over N atomic layers. This yields N — 1 
decoupled modes, corresponding to independent 
solutions of the equation V" Jm\^{z m )\ 2 = 0, 
where m numerates layers [7J. All these modes 
are, indeed, decoupled from electron system pro- 
vided that one neglects the equilibrium electron 
exchange field acting on the ion spins. In this 
approximation, we find from Eq. (JT]) that the 
modes have no dispersion and their frequencies 
coincide and are equal to f2j. In fact, weak in- 
teraction with the electrons gives rise to a small 
splitting of the ion Zeeman energies, which be- 
come dependent on the atomic layer number to: 
ftfljm = hilj + an e s% q \ 1 i l {z m )\ 2 . Taking into ac- 
count this splitting, one finds that in a symmetric 
quantum well, which we consider here, the de- 
coupled modes with antisymmetric distribution 



oc S„ 



still obey 



of ion spins J„ 
the condition J = 0, thus having no dispersion. 
For too— th mode, the ion spin precession fre- 
quency is equal to flj mo . As for the modes with a 
symmetric distribution, they become weakly cou- 
pled to the electron collective mode. One can 
show, however, that the corresponding dispersion 
is very weak provided that n e /nj <C 1. Symmet- 
ric modes also become weakly coupled to the sin- 
gle electron excitations and, consequently, slowly 
decay in the region of the Stoner continuum with 
the characteristic rate S 2 VF\/k 2 — k 2 n /tt e H,jN. 
To conclude, we have developed a theory of 



the spin waves in the 2D diluted magnetic semi- 
conductors. We have described analytically two 
collective modes corresponding to the coupled ex- 
citations of the electron and ion spins, and a large 
number of decoupled excitations of the ion spins. 
Our main finding is the tunable anticrossing of 
the collective modes. We have also predicted a 
field-induced change in a sign of the group veloc- 
ity of the ion mode and have calculated the decay 
of the waves in the Stoner continuum. 
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